Advanced ensemble learning strategy based semi-supervised soft sensing method

ABSTRACT

The present disclosure provides a novel advanced ensemble learning strategy for soft sensor development with semi-supervised model. The main target of the soft sensor is to improve the prediction performance with a limited number of labeled data samples, under the ensemble learning framework. Firstly, in order to improve the prediction accuracy of sub-models for ensemble modeling, a novel sample selection mechanism is established to select the most significantly estimated data samples. Secondly, the Bagging method is employed to both of the labeled and selected data-set, and the two different kinds of datasets are matched based on the Dissimilarity (DISSIM) algorithm. As a result, the proposed method guarantees the diversity and accuracy of the sub-models which are two significant issues of the ensemble learning. In this work, the soft sensor is constructed upon the Gaussian Process Regression (GPR) model.

TECHNICAL FIELD

The present disclosure relates to industrial modeling and soft sensing. More particularly, the present disclosure relates to the advanced ensemble learning strategy based semi-supervised soft sensing method.

BACKGROUND

Nowadays, with the growing demand of the high-quality chemical/biological products, a soft sensor has become more and more important to obtain reliable and real-time predictions of the key quality variables in industrial processes. However, it is hard to install hardware sensors to measure those key quality variables, because of the badly circumstance in the industrial environment.

In order to solve the issue that it is hard to measure key quality variables via hardware sensors, soft sensing technology has been paid more and more attentions. With the wide installation of the distributed control system (DCS) which is capable to record a large amount of process data in industrial processes, a variety of data-based soft sensors have been developed. To date, the widely used soft sensing methods include: principal component regression (PCR), partial least square (PLS), artificial neural network (ANN), support vector machine (SVM), Gaussian process regression (GPR), etc.

Nevertheless, the data samples with the output variables are difficult to be acquired, since analyzing quality variables is time-consuming and involves human resource and expensive instruments in the laboratory. As a result, we may need to corporate the soft sensor modeling problem under a very limited number of data samples for both process and quality variables, and a large amount of data samples for only process variables.

SUMMARY

In order to deal with the unbalanced numbers of the training data, this invention is motivated to build an ensemble learning based semi-supervised regression model which incorporates both labeled and unlabeled data samples.

Firstly, in order to improve the diversity performance, not only the Bagging algorithm is implemented upon both of the labeled and estimated dataset to generate sub-models, but also the way to match the labeled sub-datasets and estimated sub-datasets is designed to improve the diversity performance.

The Dissimilarity (DISSIM) algorithm is employed to measure the dissimilarity index of each labeled sub-dataset and estimated sub-dataset pair, and the most dissimilar one can be matched regularly.

Another important issue is how to guarantee the accuracy of the sub-models, which means that significant estimated data samples should be selected carefully into the next ensemble learning.

Therefore, we define an evaluation index to evaluate the important degree of each estimated data samples based on the influence that each estimated data sample worked on the labeled dataset. In this way, the estimated data samples can be selected according to the evaluation index, which means that the significant estimated data samples with additional information will be selected into the selected dataset repairing for the next ensemble step.

In summary, the missing quality variables are estimated through the K-nearest neighborhood method. Then, the important estimated samples are picked out, while the useless samples are discarded, according to the defined evaluation index. The selected samples and labeled samples are divided into several sub data sets by the Bagging algorithm. These two kinds of the sub data sets are combined based on the DISSIM algorithm. Finally, construct sub models for each combined sub data set to make a prediction for the query sample.

This invention is implemented through the following steps:

Step 1: Collect training data set which can be divided as two counterparts: 1) the labeled training data set D^(L)={X^(L)∈R^(n) ^(l) ^(×m), y^(L)∈R^(n) ^(l) ^(×1)}, 2) the unlabeled training data set D^(U)={X^(U)∈R^(n) ^(u) ^(×m),y^(U)∈R^(n) ^(u) ^(×1)};

Then use the K-nearest neighborhood method to estimate the missing quality variable in the unlabeled training data set, and denote the estimated and selected data set as D^(E)={X^(U)∈R^(n) ^(l) ^(×m), y^(E)∈R^(n) ^(l) ^(×1)} and D^(S)={X^(S)∈R^(n) ^(s) ^(×m),y^(S)∈R^(n) ^(s) ^(×1)}, respectively.

where the superscript L, U, E and S, present the labeled, unlabeled, estimated and selected data samples, respectively; and {X} and {y} are the input (process) and output (quality) variables, respectively; n_(l), n_(u), and n_(s) are the number of the labeled, unlabeled, and selected data samples, respectively; while m is the number of the process variables;

Step 2: Form the selected dataset D^(S)={X^(S)∈R^(n) ^(s) ^(×m), y^(S)∈R^(n) ^(s) ^(×1)} according to the labeled training data set and the unlabeled training data set;

where the superscript S means the selected data samples, and n_(s) is the number of the selected data samples;

Step 3: Employ Bagging algorithm to D^(L) and D^(S) respectively, then A labeled sub-datasets {D_(i) ^(L)}_(i=1) ^(A) and A selected sub-datasets {D_(j) ^(S)}_(j=1) ^(A) have been obtained; A is the number of the sub datasets;

Step 4: Match all the sub-datasets, then A matched sub-datasets have been obtained;

Step 5: When a query sample comes, predict its quality variable with the matched sub-datasets via the Gaussian Process Regression (GPR) model;

Optionally, between the step 1 and step 2, said method comprises the step:

Step 1.1: Use Eqs. (1) to (3) to calculate the estimated values of the missing output variables to form the estimated data set; said Eqs. (1) to (3) are termed as self-training algorithm;

$\begin{matrix} {{S\left( {i,j} \right)} = {{\gamma\mspace{11mu}{\exp\left( {- {{x_{i} - x_{j}}}} \right)}} + {\left( {1 - \gamma} \right)\cos\mspace{11mu}\left\langle {x_{i},x_{j}} \right\rangle}}} & (1) \\ {{{RS}\left( {i, \cdot} \right)} = {{descend}\left\{ {S\left( {i, \cdot} \right)} \right\}}} & (2) \\ {{\hat{y}}_{i} = \frac{\sum\limits_{k = 1}^{K}{{{RS}\left( {i, \cdot} \right)}_{k}y_{{{RS}{({i, \cdot})}}_{k}}^{L}}}{\sum\limits_{k = 1}^{K}{{RS}\left( {i, \cdot} \right)}_{k}}} & (3) \end{matrix}$

where ∥·∥ denotes the Euclidean distance, and γ is the parameter for the computing the similarity, and RS(i,·)_(k) represents the k-th element in the rearranged similarity vector RS(i,·), S(i,·) is the similarity vector whose elements present the similarity between the x_(i) and all the labeled samples, RS(i,·) is a vector that rearranged in descend for the similarity vector S(i,·) and y_(RS(i,·)) _(k) ^(L) is the quality variable which has the similarity RS(i,·)_(k) with the unlabeled data sample x_(i);

Optionally, said step 2 comprises the steps:

Step 2.1: Construct temporary datasets which can be denoted as D_(temp)(i)={[x^(L);x_(i) ^(U)],[y^(L);ŷ_(i) ^(E)]}. Therefore, we have n_(u) temporary datasets and they are needed for calculating the evaluation index Q; The definition of the evaluation index Q is given as follow: Q(i)=Q(x _(i) ^(U) ,ŷ _(i) ^(E))=Σ_(x) _(j) _(∈X) _(L) ({tilde over (y)} _(ji) ^(L) −y _(j) ^(L))²  (4)

where y_(j) ^(L) is the corresponding quality variable of the j-th data sample x_(j) in the labeled dataset X^(L)∈R^(n) ^(l) ^(×m), and {tilde over (y)}_(ji) ^(L) the prediction of x_(j) in the labeled dataset X^(L)∈R^(n) ^(l) ^(×m) specified by the GPR model with the i-th temporary dataset D_(temp)(i)={[x^(L);x_(i) ^(U)],[y^(L);ŷ_(i) ^(E)]};

Step 2.2: Re-arranged the index Q in a descend order; RQ=descend(Q)=RQ ₁ > . . . >RQ _(i) > . . . >RQ _(n) _(u)   (5)

where RQ_(i) represents the i-th largest value of the element in Q;

Step 2.3: Select the first n_(s) data samples with the largest Q values to form the selected dataset D^(S)={X^(S)∈R^(n) ^(s) ^(×m),y^(S)∈R^(n) ^(s) ^(×1)}.

Optionally, said step 3 comprises the following step:

Step 3.1: Calculate the dissimilarity index matrix d whose ij-th element d(i,j) represents the dissimilarity between D_(i) ^(E) and D_(j) ^(S);

In Bagging algorithm, the m samples of the l-th sub-dataset are all randomly drawn from the original dataset D, which can be represented as D_(l)={(x_(j),y_(j))}_(j=1) ^(m).

Optionally, said step 4 comprises the following steps:

Step 4.1: Search the largest value d(i_(max)) of the matrix d and match the i-th_(max) labeled sub-dataset and the j-th_(max) selected sub-dataset;

Step 4.2: Update dissimilarity index matrix d by removing the row d(i_(max)( )) and the column d(·,j_(max)( )) from the original matrix;

Step 4.3: where matrix d is the dissimilarity index matrix whose ij-th represents the dissimilarity between the i-th labeled sub-dataset and the j-th selected sub-dataset, based on the dissimilarity (DISSIM) algorithm.

Optionally, the procedure of the DISSIM algorithm of said step 4.3 is given as follows:

Consider the following two datasets, X₁∈R^(n) ¹ ^(×m) and X₂∈R^(n) ² ^(×m), where n_(i) is the number of samples in dataset X_(i) and m is the number of variables. The covariance matrices are given by

$\begin{matrix} {R_{i} = {\frac{1}{N_{i}}X_{i}^{T}X_{i}}} & (6) \end{matrix}$

The covariance matrix of the mixture of two data sets is given by

$\begin{matrix} {R = {{{\frac{N_{1}}{N_{1} + N_{2}}\begin{bmatrix} X_{1} \\ X_{2} \end{bmatrix}}^{T}\begin{bmatrix} X_{1} \\ X_{2} \end{bmatrix}} = {{\frac{N_{1}}{N_{1} + N_{2}}R_{1}} + {\frac{N_{2}}{N_{1} + N_{2}}R_{2}}}}} & (7) \end{matrix}$

Then R can be diagonalized by an orthogonal matrix P₀ using eigenvalue decomposition P ₀ ^(T) RP ₀=Λ  (8)

where Λ is a diagonal matrix whose diagonal elements are eigenvalues of R.

Then the original data X_(i) are transformed into Y_(i).

$\begin{matrix} {Y_{i} = {{\sqrt{\frac{N_{i}}{N_{1} + N_{2}}}X_{i}P_{0}\Lambda^{- \frac{1}{2}}} = {\sqrt{\frac{N_{i}}{N_{1} + N_{2}}}X_{i}P}}} & (9) \end{matrix}$

where P is transformation matrix defined as

$P = {P_{0}\Lambda^{- \frac{1}{2}}}$

The covariance matrices of the transformed data are given as

$\begin{matrix} {S_{i} = {{\frac{1}{N_{i}}Y_{i}^{T}Y_{i}} = {{\frac{N_{i}}{N_{1} + N_{2}}P^{T}\frac{X_{i}^{T}X_{i}}{N_{i}}P} = {\frac{N_{i}}{N_{1} + N_{2}}P^{T}R_{i}P}}}} & (10) \end{matrix}$

Accordingly, the index D is defined for evaluating the dissimilarity of data sets

$\begin{matrix} {D = {{{diss}\left( {X_{1},X_{2}} \right)} = {\frac{4}{m}{\sum\limits_{j = 1}^{m}\left( {\lambda^{j} - 0.5} \right)^{2}}}}} & (11) \end{matrix}$

where λ^(j) denotes the eigenvalues of the covariance matrix S_(i);

In this work, only on those eigenvalues whose values are far away from 0.5 are focused on, and the dissimilarity index {tilde over (D)} is re-defined as follow

$\begin{matrix} {\overset{\sim}{D} = {{{dis}{s\left( {X_{1},X_{2}} \right)}} = {\frac{4}{\overset{˜}{m}}{\sum\limits_{j = 1}^{\overset{\sim}{m}}\left( {{\overset{˜}{\lambda}}^{j} - {0.5}} \right)^{2}}}}} & (12) \end{matrix}$

where {tilde over (λ)}^(j) is the j-th element of the set {tilde over (λ)}^(j)={descend(|λ^(j)−0.5|)}, and {tilde over (m)} is the number of selected eigenvalues; Unlike the conventional index D measures the whole dissimilar information of the datasets, the new index {tilde over (D)} relies on the most dissimilar information between two datasets, thus {tilde over (D)} is more sensitive to different datasets.

Optionally, said step 5 comprises the steps:

Step 5.1: Make the final prediction of the quality variable by combining the sub prediction results, whose calculating method is given as follow:

$\begin{matrix} {{\hat{\gamma}}_{final} = {\sum\limits_{i = 1}^{A}{\frac{1}{A}{{\hat{\gamma}}_{i}.}}}} & (13) \end{matrix}$

Optionally, the Gaussian Process Regression model of said step 5 is described as follows:

Consider a training dataset X∈R^(n×m) and y∈R^(n×1), where {x_(i)∈R^(m)}_(i=1) ^(n) is the input data with m dimensions, and {y_(i)∈R}_(i=1) ^(n) is the output data; The aim of a regression model is to find a functional relationship y=ƒ(x|θ)+ε between input and output; Particularly, the Gaussian process regression model assumes that the regression function with zero-mean Gaussian prior distribution; The description of the function is given as follow y={ƒ(x ₁),ƒ(x ₂), . . . ,ƒ(x)}˜GP(0,K)  (14)

where K is the covariance matrix whose ij-th element is k_(i) _(j) =k(x_(i),x_(j)), and in this paper, the kernel function k(·) is the commonly used squared-exponential covariance function, which is defined as follow

$\begin{matrix} {{k\left( {x_{i},x_{j}} \right)} = {{\sigma_{f}^{2}\exp^{\{\frac{{({x_{i} - x_{j}})}^{T}{({x_{i} - x_{j}})}}{{- 2}l^{- 2}}\}}} + {\delta_{ij}\sigma_{n}^{2}}}} & (15) \end{matrix}$

where δ_(ij)=1 only when i=j, otherwise δ_(ij)=0; l is length-scale; and σ_(f) ² and σ_(n) ² are terms as signal and noise variance, respectively; An appropriate hyper-parameter set Θ=[σ_(n) ²,σ_(f) ²,l] is vital for a GPR model to make an accurate prediction about the output variable; To determine the optimal value of the hyper-parameter set, the following likelihood function should be maximized in the training process

$\begin{matrix} {L = {{{- \frac{n}{2}}{\log\left( {2\pi} \right)}} - {\frac{1}{2}{\log\left( {K} \right)}} - {\frac{1}{2}y^{T}K^{- 1}y}}} & (16) \end{matrix}$

Conjugate gradient method is used to obtain the closed solution for the above maximizing problem; An alternative way to determine the optimal value is to employ the sampling methods such as Markov chain Monte Carlo (MCMC) method, which draws samples to approximate the poster distribution of the hyper-parameters;

Once we get the optimal hyper-parameter set, the GPR model is available to make a prediction about the distribution of y_(new) corresponding for x_(new); When a query data sample x_(new) comes, according to the property of the multivariate Gaussian distribution, the posterior distribution of the output, namely (y_(new)|X,y,x_(new))˜N(μ(y_(new)),σ²(y_(new))) can be identified, where μ(y_(new)) and σ²(y_(new)) is the mean and variance of the poster distribution, respectively;

They can be calculated as follows μ(y _(new))=k ^(T)(x _(new) ,X)K ⁻¹ y  (17) σ²(y _(new))=k(x _(new) ,x _(new))−k ^(T)(x _(new) ,X)K ⁻¹ k(x _(new) ,X)  (18)

where k(x_(new),X)=[k(x_(new),x₁),k(x_(new),x₂), . . . , k(x_(new),x_(n))]^(T); Finally, take the expectation μ(y_(new)) of the posterior distribution as the prediction result about the output variable of the query data sample x_(new).

Optionally, the above advanced ensemble learning strategy based semi-supervised soft sensing method is characterized by improving the prediction performance with a limited number of labeled data samples, under the ensemble learning framework, and the improved prediction performance can be obtained.

BRIEF DESCRIPTION OF FIGURES

FIG. 1 shows the modeling steps of the advanced ensemble learning strategy based semi-supervised soft sensing method in example 1;

FIG. 2 shows the block diagram of the debutanizer column process;

FIG. 3 shows the modeling steps of the advanced ensemble learning strategy based semi-supervised soft sensing method;

FIG. 4 shows the prediction errors under different modeling strategies;

FIG. 5 shows the soft sensor based control strategy of the debutanizer column process.

DETAILED DESCRIPTION Example 1

Based on FIG. 1, the detailed illustration of the advanced ensemble learning strategy based semi-supervised soft sensing method is given.

Step 1: Collect training data set which can be divided as two counterparts: 1) the labeled training data set D^(L)={X^(L)∈R^(n) ^(l) ^(×m),y^(L)∈R^(n) ^(l) ^(×1)}, 2) the unlabeled training data set D^(U)={X^(U)∈R^(n) ^(u) ^(×m),y^(U)∈R^(n) ^(u) ^(×1)};

Then use self-training method to estimate the missing quality variable in the unlabeled training data set, and denote the estimated and selected data set as D^(E)={X^(U)∈R^(n) ^(l) ^(×m),y^(E)∈R^(n) ^(l) ^(×1)} and D^(S)={X^(S)∈R^(n) ^(s) ^(×m),y_(S)∈R^(n) ^(s) ^(×1)}, respectively.

where the superscript L, U, E and S, present the labeled, unlabeled, estimated and selected data samples, respectively; and {X} and {y} are the input (process) and output (quality) variables, respectively; n_(l),n_(u) and n_(s) are the number of the labeled, unlabeled, and selected data samples, respectively; while m is the number of the process variables;

Step 2: Form the selected dataset D^(S)={X^(S)∈R^(n) ^(s) ^(×m),y^(s)∈R^(n) ^(s) ^(×1)} according to the labeled training data set and the unlabeled training data set;

where the superscript S means the selected data samples, and n_(s) is the number of the selected data samples;

Step 3: Employ Bagging algorithm to D^(L) and D^(S) respectively, then A labeled sub-datasets {D_(i) ^(L)}_(i=1) ^(A) and A selected sub-datasets {D_(j) ^(S)}_(j=1) ^(A) have been obtained; A is the number of the sub datasets;

Step 4: Match all the sub-datasets, then A matched sub-datasets have been obtained;

Step 5: When a query sample comes, predict its quality variable with the matched sub-datasets via the Gaussian Process Regression (GPR) model.

Example 2

Based on FIG. 1, the detailed illustration of the advanced ensemble learning strategy based semi-supervised soft sensing method is given.

The debutanizer column process is utilized, and the soft sensors are designed for the butane content. The block diagram of the debutanizer column process is demonstrated in FIG. 2. In this process, the content of butane in the bottom of the column is the primary variable which should be predicted in real-time. Seven hardware sensors which are marked with grey circles in FIG. 2, are installed to obtain secondary variables, the explanations of which are presented in Table 1.

TABLE 1 Secondary Secondary variable Explanation variable Explanation x₁ Top x₅ 6th tray temperature temperature x₂ Top x₆ Bottom pressure temperature I x₃ Reflux x₇ Bottom flow temperature II x₄ Flow to next process

Step 1: Collect training data set which can be divided as two counterparts:

1) the labeled training data set D^(L)={X^(L)∈R^(n) ^(l) ^(×m),y^(L)∈R^(n) ^(l) ^(×1)};

2) the unlabeled training data set D^(U)={X^(U)∈R^(n) ^(u) ^(×m),y^(U)∈R^(n) ^(u) ^(×1)}.

Then use self-training method to estimate the missing quality variable in the unlabeled training data set, and denote the estimated and selected data set as D^(E)={X^(U)∈R^(n) ^(l) ^(×m),y^(E)∈R^(n) ^(l) ^(×1)} and D^(S)={X^(S)∈R^(n) ^(s) ^(×m),y^(S)∈R^(n) ^(s) ^(×1)}, respectively.

where the superscript L, U, E and S present the labeled, unlabeled, estimated and selected data samples, respectively; and {X} and {y} are the input (process) and output (quality) variables, respectively; n_(l) and n_(u) are the number of the labeled and unlabeled data samples, respectively; while m is the number of the process variables.

Then use Eqs. (1) to (3) to calculate the estimated values of the missing output variables to form the estimated data set, Eqs. (1) to (3) are termed as self-training algorithm.

$\begin{matrix} {{S\left( {i,j} \right)} = {{{\gamma exp}\left( {- {{x_{i} - x_{j}}}} \right)} + {\left( {1 - \gamma} \right)\cos\left\langle {x_{i},x_{j}} \right\rangle}}} & (1) \\ {{{RS}\left( {i, \cdot} \right)} = {{descend}\left\{ {S\left( {i, \cdot} \right)} \right\}}} & (2) \\ {{\hat{y}}_{i} = \frac{\sum\limits_{k = 1}^{K}{{{RS}\left( {i, \cdot} \right)}_{k}y_{R{S{({i, \cdot})}}_{k}}^{L}}}{\sum\limits_{k = 1}^{K}{R{S\left( {i, \cdot} \right)}_{k}}}} & (3) \end{matrix}$

where ∥·∥ denotes the Euclidean distance, and γ is the parameter for the computing the similarity, and RS(i,·)_(k) represents the k-th element in the rearranged similarity vector RS(i,·), and y_(RS(i,·)) _(k) ^(L) is the quality variable which has the similarity RS(i,·)_(k) with the unlabeled data sample x_(i).

Step 2: Then construct temporary datasets which can be denoted as D_(temp)(i)={[x^(L);x_(i) ^(U)],[y^(L);ŷ_(i) ^(E)]}. Therefore, we have n_(u) temporary datasets and they are needed for calculating the evaluation index Q. The definition of the evaluation index Q is given as follow Q(i)=Q(x _(i) ^(U) ,ŷ _(i) ^(E))=Σ_(x) _(j) _(∈X) _(L) ({tilde over (y)} _(j) ^(L) −y _(j) ^(L))²  (4)

where y_(j) ^(L) is the corresponding quality variable of the j-th data sample x_(j) in the labeled dataset X^(L)∈R^(n) ^(l) ^(×m), and {tilde over (y)}_(j) ^(L) is the prediction of x_(j) specified by the GPR model with the i-th temporary dataset D_(temp)(i) {[x^(L);x_(i) ^(U)],[y^(L);ŷ_(i) ^(E)]}.

Re-arranged the index Q in a descend order RQ=descend(Q)=RQ ₁ > . . . >RQ _(i) > . . . >RQ _(n) _(u)   (5)

where RQ_(i) represents the i-th largest value of the element in Q.

Select the first n_(s) data samples with the largest Q values to form the selected dataset D^(S)={X^(S)∈R^(n) ^(s) ^(×m),y^(s)∈R^(n) ^(s) ^(×1)}.

Step 3: Employ Bagging algorithm to D^(L) and D^(S) respectively, then A labeled sub-datasets {D_(i) ^(L)}_(i=1) ^(A) and A selected sub-datasets {D_(j) ^(S)}_(j=1) ^(A) have been obtained. Calculate the dissimilarity index matrix d whose ij-th element d(i,j) represents the dissimilarity between D_(i) ^(E) and D_(j) ^(S).

In Bagging algorithm, the m samples of the l-th sub-dataset are all randomly drawn from the original dataset D, which can be represented as D_(l)={(x_(j),y_(j))}_(j=1) ^(m).

Step 4: Repeat this step until all the sub-datasets are matched. Search the largest value d(i_(max)) of the matrix d and match the i-th_(max) labeled sub-dataset and the j-th_(max) selected sub-dataset. Then update dissimilarity index matrix d by removing the row d(i_(max)( )) and the column d(·,j_(max)( )) from the original matrix.

where matrix d is the dissimilarity index matrix whose ij-th represents the dissimilarity between the i-th labeled sub-dataset and the j-th selected sub-dataset, based on the dissimilarity (DISSIM) algorithm.

The procedure of the DISSIM algorithm is given as follows:

Consider the following two datasets, X₁∈R^(n) ¹ ^(×m) and X₂∈R^(n) ² ^(×m), where n_(i) is the number of samples in dataset X_(i) and m is the number of variables. The covariance matrices are given by

$\begin{matrix} {R_{i} = {\frac{1}{N_{i}}X_{i}^{T}X_{i}}} & (6) \end{matrix}$

The covariance matrix of the mixture of two data sets is given by

$\begin{matrix} {R = {{{\frac{N_{1}}{N_{1} + N_{2}}\begin{bmatrix} X_{1} \\ X_{2} \end{bmatrix}}^{T}\begin{bmatrix} X_{1} \\ X_{2} \end{bmatrix}} = {{\frac{N_{1}}{N_{1} + N_{2}}R_{1}} + {\frac{N_{2}}{N_{1} + N_{2}}R_{2}}}}} & (7) \end{matrix}$

Then R can be diagonalized by an orthogonal matrix P₀ using eigenvalue decomposition P ₀ ^(T) RP ₀=Λ  (8)

where Λ is a diagonal matrix whose diagonal elements are eigenvalues of R.

Then the original data X_(i) are transformed into Y_(i)

$\begin{matrix} {Y_{i} = {{\sqrt{\frac{N_{i}}{N_{1} + N_{2}}}X_{i}P_{0}\Lambda^{- \frac{1}{2}}} = {\sqrt{\frac{N_{i}}{N_{1} + N_{2}}}X_{i}P}}} & (9) \end{matrix}$

where P is transformation matrix defined as

$P = {P_{0}\Lambda^{- \frac{1}{2}}}$

The covariance matrices of the transformed data are given as

$\begin{matrix} {S_{i} = {{\frac{1}{N_{i}}Y_{i}^{T}Y_{i}} = {{\frac{N_{i}}{N_{1} + N_{2}}P^{T}\frac{X_{i}^{T}X_{i}}{N_{i}}P} = {\frac{N_{i}}{N_{1} + N_{2}}P^{T}R_{i}P}}}} & (10) \end{matrix}$

Accordingly, the index D is defined for evaluating the dissimilarity of data sets

$\begin{matrix} {D = {{dis{s\left( {X_{1},X_{2}} \right)}} = {\frac{4}{m}{\sum\limits_{j = 1}^{m}\left( {\lambda^{j} - {0.5}} \right)^{2}}}}} & (11) \end{matrix}$

where λ^(j) denotes the eigenvalues of the covariance matrix S_(i).

In this work, only on those eigenvalues whose value is far away from 0.5 is focused on, and the dissimilarity index {tilde over (D)} is re-defined as follow

$\begin{matrix} {\overset{\sim}{D} = {{{dis}{s\left( {X_{1},X_{2}} \right)}} = {\frac{4}{\overset{˜}{m}}{\sum\limits_{j = 1}^{\overset{\sim}{m}}\left( {{\overset{˜}{\lambda}}^{j} - {0.5}} \right)^{2}}}}} & (12) \end{matrix}$

where {tilde over (λ)}^(j) is the j-th element of the set {tilde over (λ)}^(j)={descend(|λ^(j)−0.5|)}, and {tilde over (m)} is the number of selected eigenvalues. Unlike the conventional index D measures the whole dissimilar information of the datasets, the new index {tilde over (D)} relies on the most dissimilar information between two datasets, thus {tilde over (D)} is more sensitive to differentiate datasets.

Step 5: Based on the above steps, the A sub-datasets have been obtained. When a query sample comes, predict its quality variable with the sub-datasets via the Gaussian Process Regression (GPR) model. Make the final prediction of the quality variable by combining the sub prediction results, whose calculating method is given as follow

$\begin{matrix} {{\hat{y}}_{final} = {\sum\limits_{i = 1}^{A}{\frac{1}{A}{\hat{y}}_{i}}}} & (13) \end{matrix}$

The Gaussian Process Regression model is described as follows:

Consider a training dataset X∈R^(n×m) and y∈R^(n×1), where {x_(i)∈R^(m)}_(i=1) ^(n) is the input data with m dimensions, and {y_(i)∈R}_(i=1) ^(n) is the output data. The aim of a regression model is to find a functional relationship y=ƒ(x|θ)+ε between input and output. Particularly, the Gaussian process regression model assumes that the regression function with zero-mean Gaussian prior distribution. The description of the function is given as follows y={ƒ(x ₁),ƒ(x ₂), . . . ,ƒ(x _(n))}˜GP(0,K)  (14)

where K is the covariance matrix whose ij-th element is k_(ij)=k(x_(i),x_(j)), and in this paper, the kernel function k(·) is the commonly used squared-exponential covariance function, which is defined as follow

$\begin{matrix} {{k\left( {x_{i},x_{j}} \right)} = {{\sigma_{f}^{2}\exp^{\{\frac{{({x_{i} - x_{j}})}^{T}{({x_{i} - x_{j}})}}{{- 2}l^{- 2}}\}}} + {\delta_{ij}\sigma_{n}^{2}}}} & (15) \end{matrix}$

where δ_(ij)=1 only when i=j, otherwise δ_(ij)=0; l is length-scale; and σ_(f) ² and σ_(n) ² are terms as signal and noise variance, respectively. An appropriate hyper-parameter set Θ=[σ_(n) ²,σ_(f) ²,l] is vital for a GPR model to make an accurate prediction about the output variable. To determine the optimal value of the hyper-parameter set, the following likelihood function should be maximized in the training process

$\begin{matrix} {L = {{{- \frac{n}{2}}{\log\left( {2\pi} \right)}} - {\frac{1}{2}{\log\left( {K} \right)}} - {\frac{1}{2}y^{T}K^{- 1}y}}} & (16) \end{matrix}$

Conjugate gradient method is used to obtain the closed solution for the above maximizing problem. An alternative way to determine the optimal value is to employ the sampling methods such as Markov chain Monte Carlo (MCMC) method, which draws samples to approximate the poster distribution of the hyper-parameters.

Once we get the optimal hyper-parameter set, the GPR model is available to make a prediction about the distribution of y_(new) corresponding for x_(new). When a query data sample x_(new) comes, according to the property of the multivariate Gaussian distribution, the posterior distribution of the output, namely (y_(new)|X,y,x_(new))˜N(μ(y_(new)),σ²(y_(new))) can be identified, where μ*y_(new)) and σ²(y_(new)) is the mean and variance of the poster distribution, respectively, which can be calculated as follows μ(y _(new))=k ^(T)(x _(new) ,X)K ⁻¹ y  (17) σ²(y _(new))=k(x _(new) ,x _(new))−k ^(T)(x _(new) ,X)K ⁻¹ k(x _(new) ,X)  (18)

where k(x_(new),X)=[k(x_(new),x₁),k(x_(new),x₂), . . . , k(x_(new),x_(n))]^(T). Finally, take the expectation μ(y_(new)) of the posterior distribution as the prediction result about the output variable of the query data sample x_(new).

FIG. 4 shows the prediction errors under different modeling strategies, it can be found that the advanced ensemble learning strategy based semi-supervised soft sensing method is superior to other strategies.

The main operating variable affecting the product quality is the reflux flow x₃. Thus, a soft sensor based control strategy is designed to ensure the butane concentration of the product by manipulating the reflux rate. The close loop control scheme is shown in FIG. 5.

The signals of secondary variables are smoothed before sending into the soft sensor model. Furthermore, in order to make the control system ensure both the product quality and the material balance, the “expert control module” is applied in the inference control strategy. If the butane concentration is qualified, the reflux rate shall be adjusted according to x₂ and x₄; otherwise, the reflux rate shall be adjusted according to the deviation between the butane estimated value and the set value in the soft instrument. Finally, set the upper and lower bounds of the valve opening according to the control requirements of the distillation column.

While various aspects and embodiments have been disclosed herein, other aspects and embodiments will be apparent to those skilled in the art. The various aspects and embodiments disclosed herein are for purposes of illustration and are not intended to be limiting, with the true scope and spirit being indicated by the following claims. 

What is claimed is:
 1. A method comprising: obtaining a labeled training set of vectors of process parameters of a debutanizer column process and labels respectively associated with the vectors of the labeled training set; obtaining an unlabeled training set of vectors of the processor parameters of the debutanizer column process; obtaining an estimated data set by labeling the vectors of the unlabeled training set utilizing a self-training method; selecting a subset of the estimated data set as a selected data set, based on the labeled training set and the unlabeled training set; obtaining subsets of the labeled training set and obtaining subsets of the selected data set, utilizing a bagging algorithm; pairing the subsets of the labeled training set and the subsets of the selected data set; calculating a butane content at a bottom of a debutanizer column used in the debutanizer column process, from a vector of the processor parameters of the debutanizer column process based on the paired subsets of the labeled training set and subsets of the selected data set; and adjusting the process parameters of the debutanizer column process based on the butane content; wherein the process parameters are selected from the group consisting of a temperature at a top of the debutanizer column, a pressure at the top of the debutanizer column, a reflux flow rate in the debutanizer column, a flow rate to a downstream process, tray temperatures and combinations thereof.
 2. The method of claim 1, wherein adjusting the process parameters of the debutanizer column process based on the butane content comprises maximizing the butane content.
 3. The method of claim 1, wherein selecting the subset of the estimated data set comprises: for each member of the estimated data set, calculating a sum of squares of differences between labels associated with the vectors of the labeled training set and a label of that member; and selecting the subset of the estimated data set in a matter that the sum calculated for any member of the subset of the estimated data set is not smaller than the sum calculated for any member of the estimated data set that is not of the subset of the estimated data set.
 4. The method of claim 1, wherein obtaining the subsets of the labeled training set and obtaining the subsets of the selected data set comprise calculating a dissimilarity index for each pair of a member of the labeled training set and a member of the selected data set.
 5. The method of claim 4, wherein pairing the subsets of the labeled training set and the subsets of the selected data set is based on the dissimilarity index. 